Loss of superfluidity by fermions in the boson Hubbard model on an optical lattice 
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The experimentally observed loss of superfluidity by introducing fermions to the boson Hubbard 
system on an optical lattice is explained. We show that the virtual transitions of the bosons to the 
higher Bloch bands, coupled with the contact boson-fermion interactions of either sign, result in 
an effective increase of the boson on-site repulsion. If this renormalization of the on-site potential 
is dominant over the fermion screening of the boson interactions, the Mott insulating lobes of the 
Bose-Hubbard phase diagram will be enhanced for either sign of the boson-fermion interactions. We 
discuss implications for cold atom experiments where the expansion of the Mott lobes by fermions 
has been conclusively established. 
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Introduction: The intrinsic effect of introducing a gas 
of fermions on the famous superfluid-insulator (SI) phase 
diagram [l[ of the boson Hubbard model is still theoret- 
ically unresolved. Experimentally, the situation is quite 
unambiguous 0, 0, 0|: the addition of a degenerate gas 
of spin-polarized fermions to the bosonic superfluid con- 
densate reduces the superfluid coherence, irrespective of 
the sign of the interaction between the bosons and the 
fermions. This indicates that the condensate, in the pres- 
ence of interactions with the fermions, gives way to the 
Mott insulating phase at values of the effective t/U larger 
than that without the fermions, where t is the nearest 
neighbor hopping amplitude and U is the on-site inter- 
action of the constituent bosons. The areas of the Mott 
insulating lobes in the boson Hubbard phase diagram are 
thus enhanced by the fermions. The earlier experiments 
0, Q observed the loss of superfluid coherence for fixed 
attractive boson-fermion interactions, Ubfj which were 
larger in magnitude than the boson on-site repulsion it- 
self. Recently this finding has also been confirmed for 
both attractive and repulsive interspecies interactions in 
a range of values for |!7bf| both smaller and larger than 
U [3| ■ While the problem of the Bose- Fermi mixture is 
theoretically interestingfor the ens uing complexity in the 
phase diagram [1, d, 0, H, S 13, U, 12, 13, EH, the main 
experimental observation - a suppression of the super- 
fluid coherence by adding fermions - is contrary to many 
recent theoretical studies [3, 12, 14 1; the only notable 



exception is the numerical work in Ref. [la ]. 

In the theoretical studies of the SI phase diagram of the 
Bose-Fermi mixture SBEEEl, the central argu- 
ment involves the screening of the boson on-site repulsive 
potential by the fermions. The mobile fermions mediate a 
spatially and temporally non-local attractive interaction 
among the bosons, which can be shown to screen and re- 
duce the repulsive interaction U of the boson Hubbard 
model. Since the on-site repulsion tends to localize the 
bosons on the individual lattice sites (favoring the Mott 



insulating phases for any integer filling factor), any re- 
duction of this interaction, such as that provided by the 
fermionic screening, should expand the area occupied by 
the superfluid phase in the phase diagram @, [H, 14|. As 
pointed out before, this simple intuitive picture, made 
rigorous recently both numerically 12j and analytically 
[lj] on the single band boson Hubbard model, cannot 
explain the available experiments. Here we show that 
there is an additional effect of the fermion contact inter- 
actions of either sign, mediated by virtual transitions of 
the bosons to the higher Bloch bands, which leads to an 
effective increase of the boson on-site interactions. There 
is some numerical evidence of this effect [13] for the case 
of attractive interspecies interactions only. In this Letter, 
we treat the above two disparate effects within a unified 
analytical framework. We believe that this theory, which 
includes the idea of virtual transitions to the higher bo- 
son Bloch bands, provides an explanation for the loss 
of bosonic superfluid coherence by introducing fermions 
(irrespective of the sign of the interspecies interactions), 
seen in the recent cold atom experiments. Our conclu- 
sion is that multi-band effects, neg lected in the previous 
studies [!, @, 0, 0, S E3, ED, El, EE E3 , are important to 
describe the physics of Bose-Fermi mixtures. 

To get a unified description of the effects of the 
fermions, we start with the multi-band boson Hubbard 
model where the two lowest-lying boson Bloch bands are 
separated by an energy gap fl t, (X, U, where [i is the 
boson chemical potential. We assume that the fermions 
are coupled to the bosons in the two lowest-lying Bloch 
bands with contact interactions. The trivial effects of the 
presence of the higher boson bands are renormalizations 
of the bare bosonic parameters t, ji, U , which are present 
even without the fermions (and hence, will be mostly 
suppressed here). However, most importantly, we find 
that the virtual transitions of the bosons to the higher 
Bloch bands give rise to a new type of on-site repulsive 
interaction among the bosons mediated by the fermions. 
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This interaction, which is nonlinear in the boson-fermion 
scattering length obf, and remains repulsive for either 
sign of p , tends to hinder the flow of the bosons on the 
optical lattice. Thus, this interaction moves the Mott in- 
sulating transitions to shallower lattice depths for any 
boson integer filling factor. Including this new interac- 
tion and the usual fermion-mediated scree ning interac- 
tion 0, HI, 14 1 in a unified framework 14 1 , we find 
that the question of the overall shift of the phase diagram 
is a quantitative one arising from a competition between 
attractive and repulsive terms: if the fermion-induced 
higher-band renormalization of the on-site potential is 
dominant, as may be the case in the recent experiments, 
superfluid coherence will be reduced by introducing the 
fermions. 

Model: We begin with the following second quantized 
Hamiltonian describing the bosons and spin-polarized 
fermions interacting with each other through a contact 
interaction in an optical lattice: 

H = j d 3 r^{r) [fl-jf + ^$ t (r)$(r)] $(r)+ (1) 



9BF 



d d r& (r)^(r)^ (r)*(r). 



Here, $(r) and fy(r) are the field operators for the 
bosonic and the fermionic atoms, respectively, the single 
particle Hamiltonians H^ F = Tb/f + Vi a t(r), with Tb/f 
representing the kinetic energy for the bosons/fermions 
and Vi a t(r) denoting the lattice potential, V\ At {r) = 
J>2j—i Vq sin 2 (tt-^) with a as the lattice spacing. The in- 
teraction coupling constants are given by gss = ^ B > 
qbf = , where tub is the mass of a bosonic atom, 

TO rc d is the boson-fermion reduced mass and asB /bf are 
the boson-boson and boson-fermion scattering lengths, 
respectively. 

We now expand the field operators in the Wannier 
function basis, $(r) = J^i a bi, a w a (r — r»), and 9(r) = 



where the operators &i jQ and Ci, 



annihilate the bosons and the fermions at a site i in 
a band a, respectively. Substituting the above expres- 
sions for the field operators into Eq. (JTJ) , one obtains 
the full multi-band model 1|| 17 , [l8| for the Bose- Fermi 
system. Instead of treating the complexity of the full 
multi-band Hamiltonian, we consider, for simplicity, a 
two-band model for the bosons with the fermions being 
in a single band, and keep the largest band-mixing terms 
in the bosonic part of the Hamiltonian. This effective 
two-band model for the bosons, which can be justified for 
large interband energy separation f2 = \ZAErVq, where 
Er = it 2 /2m b a 2 is the recoil energy, captures the essen- 
tial physics involving the virtual transitions of the bosons 
to the higher Bloch bands. The fermion Wannier wave- 
functions u a (r — Ti) are chosen using the mean-field one- 
body Hamiltonian for the fermions, Tp + Vp(r), where 



Pb{t) = n \wi,i(r)\ 2 with ps(r) and n being the aver- 
age boson density per site and average boson number 
per site, respectfully. Thus, the shapes of these functions 
within a unit cell, which will be important later to deter- 
mine the sign of the fermion renormalization of the on- 
site bosonic potential (see the discussion after Eq. ©), 
depend on the sign of the interspecies interactions (sign 
of obf). In the rest of the paper we will study the two- 
band Bose-Fermi model defined by the Hamiltonian, 



H — Hi + H[ + Ht + Hb f + Hp 



U a ' a . , A 
s a n ha + 2 n i!a {n i<a 1) 



(2) 

2,CK>Q( / 



H t =- Y, t a [4, a b ] ,a + il-c];HBF=Y U FB[<c l ~(n l ^}i 
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Hf=Y \_ e ^i S i3 ~ f F (4 C 3+ H -' 



the effective potential Vf(v) = Vi a t(r)- 
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Here, 



Here, a = 1, 2, the energies e a = (wi iCt \HQ \wi >a ) — n, the 
matrix elements U a ' a = gBB(wi,a',Wi ia >\\wi^ a ;Wi ia >) and 
Up B = gBF/2(u>i,a', Ui\wi_ a ; Ui). The fermion energy and 
hopping are given by eo = (ui\Tp + Vp(r)\ui) — /i and 
tp = —(ui\Tp + Vp{r)\uj). Note that the piece H[ in 
Eq. ((2|), which corresponds to scattering of two bosons be- 
tween the first and the second Bloch bands, leads to band 
mixing [20] . As we show below, these band-mixing terms, 
coupled with the fermion contact interactions, renormal- 
ize the local repulsive interaction between the constituent 
bosons in the low energy subspace. 

Schrieffer- Wolff transformation: In order to reveal 
the nature of the renormalization, one needs to diagonal- 
ize the full Hamiltonian and project it on the subspace 
of the lower band only. Assuming that fl is the largest 
energy scale in the problem, we perform the Schrieffer- 
Wolff canonical transformation [19( on the Hamiltonian 
H to decouple the diagonal (boson number-conserving) 
and the non-diagonal (band-mixing) pieces of the Hamil- 
tonian to a given order in 1/51: 

H cS = e s He- s = H+[S,H} + l[S, [S, H]} + ... (3) 

To look for the right unitary transformation operator S, 
we separate the Hamiltonian in Eq. ^ into three parts: 
H Q = nj2 t ni,2, H 1 = H-H -H 2 and H 2 = H[. To 
the zeroth order in the operator is given by 

5(1) =E ^(OL^'^'- & t* & tA«^«) ( 4 ) 

i.a>a.' 

One can check that S (1 > satisfies [S^,H ] = -H 2 , 
so, to the zeroth order in I/O, the effective Hamilto- 
nian is simply the sum of Hq and H±. The next order 
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operators and can be found recursively, e.g., 
g( 2 ) = [H2, Hi]/4Q 2 . Using the canonical transformation 
denned in Eq. © with S = S« + + S( 3 \ we decou- 
ple H to the order 1 /tt 2 : H eS =H +H 1 + \ [S^ , H 2 ] + 
^[S^ 2 ' , H2]. Then, by projecting H c g on the lowest Bloch 
band, and omitting the band index, we arrive at the fol- 
lowing low-energy Hamiltonian to this order: 



i 



1 fj PU FB F (U^\\ _ 



{hi))hl 



Hp 



(5) 



The tilde indicates the renormalization of the original 
boson Hubbard parameters U, t and \x due to processes 
involving virtual transitions of the bosons to the higher 
Bloch band, e.g., U = U- (U 1 ' 2 ) 2 /ft + .... These renor- 
malizations are also present in a pure bosonic system and 
are independent of Ufb- Most importantly, however, the 
contact interactions with the fermions provide an addi- 
tional renormalization to the boson on-site potential U 
which is linear in U fb ■ 

The correction to the boson-boson repulsion due to 
the fermions corresponds to a virtual process when two 
bosons are excited from the first to the second Bloch 
band. The probability for such processes to occur {i.e., 
the fraction of time the system dwells in such a virtual 
state) is proportional to (t/ 1 ' 2 /J7) 2 ni(n J — 1). While in the 
virtual state, the interaction energy between the bosons 
and the fermions (defined by Hbf in Eq. ([2])) changes by 
2(Ug F — Ug F )npi. Thus, by combining the above two 
terms and introducing the dimensionless parameter p (to 
account for the change in the interaction energy), 

p = l-(wix, Ui\wiX,Ui}/(wiX,Ui\w it i;Ui), (6) 

one recovers the correction to the boson-boson repulsion 
given in Eq. ([5]). For attractive boson- fcrmion interac- 
tions, the fcrmionic wavefunction u(r) is peaked at the 
center. Thus, the overlap of Ui with the boson Wannicr 
function in the second Bloch band, Wi t 2, is smaller than 
its overlap with Wi i. Therefore, for negative Ubf , P is 
positive and 0(1). However, p changes sign for repulsive 
Ubf, since, in this case, the two species of atoms maxi- 
mize the distance between them (i.e., the fermion density 
is suppressed at the center of the unit cell), resulting in 
the numerator in the second term in Eq. ([5]) exceeding 
the denominator. From Eq. ([5]), notice that the sign of 
the renormalization to U is determined by sgn(p[/pB)- 
Therefore, it remains repulsive and tends to suppress the 
superfluid phase for either sign of Ubf- 

Shift of the phase diagram: To the lowest order in 
Ufb, the on-site interaction energy and the chemical 
potential of the bosons are modified as U — ► U' = 



U - pU FB (U^/Vzty n° Fi and fx -> fj,' = p, - U FB n 
with n Fi the average density of the fermions. Since the 





Fii 



fermions appear only at the quadratic order in Eq. ([5]) , in- 
tegrating them out leads to the effective imaginary-time 
action 



U' 

b*d T bi + —ni(ni — l) — (in-, 



(7) 



[211 

rP 

S ee (b*,b) = V/ dr 
i J o 

- / drHt-yJ dri / dT 2 ni(T 1 )Mi j (Ti-T2)n j (T 2 ) 
Jo fj . Jo Jo 

From the last term in Eq. ([7]), it is clear that the 
fermions induce another boson-boson interaction with 
the spatially and temporally non-local kernel (to the or- 

- U 2 



der U FB 



/A) Mij(n - 75) = U FB (Sn Fi (Ti)5n Fj (T2))/2. 



Here, Sn F i 



n F i 



and the parameter A is pro- 



portional to the inverse of the density of states at the 
Fermi level v F , A = \ jvpo? . This contribution describes 
the screening of the bosonic repulsive interactions by the 
fermions, which effectively reduces U, leadin g to the sup- 
pression of the Mott insulating phase 

BETH. 

The effect of the above two competing contributions 
on the phase diagram can be calculated analytically. We 
first need to calculate the boson on-site Green's function 
22] for the action in Eq. ([7]) at zero frequency [T3 |. 



Gi(0): 



(8) 



n + l 

8E,, 



p 



pU FB n n° Fi /U 1 ' 2 



2SE P V 



A5E„ UE 



8E, 



SE h 



pU F b{1 - n Q )n Fi (U 



25E h 



7-1,2 



A5E h 



Here 5E P and SEh are the particle and the hole excitation 
energies: SE P = U'uq — fi' and SEh = \J — U'(no — 1), 
no is the number of bosons per site minimizing the 
ground state energy. The dimensionless function R(y) = 
-^J^x 2 dx r o °° di/p^n(x, v)/v F , where II(g,fi„) is 
the three-dimensional fermion polarization function [23| . 
The dependence of R(y) on its argument is plotted in 
Fig. [TJ The mean field SI phase boundary can be ob- 
tained by solving the equation l/zt+ J^drG^T) = 0, 
where z is the lattice coordination number [23| • 

To the linear order in U F b , the correction to the phase 
boundary is influenced only by the enhancement of the 
on-site repulsion given in Eq. ([5J): 



5tM _t(a B F)-t(a B F = 0) UpBpfU 1 



U 



U 



2U V fi 



(9) 



[1 + 2(p/U) 2 - 2{fj/U)(n - 1) - n }n 
[1 + WU)} 2 



Notice that for either attractive or repulsive Ubf, the 
product of Ubf and p remains negative. Therefore, the 
shift St^ is towards shallower lattices, indicating an ex- 
pansion of the Mott-insulating lobes. The shift of the 
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FIG. 1: (Color online) (a) Main panel: Phase Boundary of 
the boson Hubbard model with and without the fermions for 
the boson density no = 5. Solid line describes the insulator- 
superfluid phase boundary without the fermions. The dashed 
line corresponds to the phase boundary with the fermions 
present in the region of the validity of the perturbation se- 
ries [H]. Here we used Ufb ~ -Ubb, U 1 ' 2 /fi w 0.3, 
A/Ubb = 40 and n° Fi = 0.75. (b) Inset: The dependence 
of the function R(y) on its argument. 



phase boundary due to the fermion-mediated screening, 
which manifests itself only in the second order in Ufb, 
is given by, 



U 2 

u FB 



u 



AU 



H)' 



n n 




+(n +l)l^-(no-l) 



As expected, the fermion-mediated screening enhances 
the area occupied by the superfluid phase. Notice that 
the small perturbation parameters used in the above cal- 
culations are an d ^f~, respectively. Therefore, they 
can be parametrically valid for both | Ub f | smaller as well 
as larger than the bare U. At the tip of the Mott lobes, 
the ratio of the two contributions for rig 1 is given by, 



5tV 



5t( 2 ) 



]p\n° F A 



\U FB \noR(^) 



U 12 

n 



(ii) 



Since the two competing contributions depend on differ- 
ent independent parameters, the sign of the phase bound- 
ary shift is a quantitative question. For example, for 
(■^-) 2 ~ U ^ B ^ , the two effects are comparable. In gen- 



eral, for 



\u FE 



< 



the superfluid state is suppressed 



a ^ n ' 

for either sign of the interspecies interaction, as shown in 
the figure Q] 

Conclusion: We show that the virtual transitions of 
bosons to the higher Bloch bands in an optical lattice, 



coupled with contact interactions with a degenerate gas 
of fermions, generate a new renormalization of the in- 
teractions in the boson Hubbard model. For either sign 
of the coupling between the fermions and the bosons, 
this renormalization enhances the boson on-site repul- 
sion, and thus favors the Mott insulating phase. If this ef- 
fect is dominant over the usual fcrmion mediated screen- 
ing, the superfluid coherence of the Bose-Hubbard system 
will be suppressed by the fermions, as has been observed 
in recent experiments [H, 0, 0] • Finally, we emphasize that 
promoting bosons to higher Bloch bands can lead to new 
types of heteronuclear interactions which are important 
for the realistic description of heteronuclear mixtures in 
cold-atom experiments. 
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